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ABSTRACT 


Numerical solutions are obtained for two-dimensional, unsteady 
potential flow about a flat plate. Two types of flow are investigated: 
small amplitude oscillation about a spanwise axis and impulsive accelera- 
tion of a nonoscillating plate. Pressure, forces and moments, and kine- 
matics of the vortex wake are calculated. 

In this work vortex sheet shedding from the trailing edge is 
included. In future efforts, vortex sheet shedding from trailing and 


leading edges are contemplated. 





NOMENCLATURE 
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1 [pee 


<< 


nM >in > 2 
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chord length 

drag per unit span 

unit vectors in (x-, y-, z-) directions respectively 
lift per unit span 

number of time increments taken 

number of evenly distributed bound vortices 
position vector 

arbitrary running variable 

time 

trailing vortex sheet 

free stream velocity 

final free stream velocity 

field velocity at ge 

vortex sheet 


absolute reference system 


vortex strength distribution 

free vortex strength 

amplitude of oscillation 

body fixed reference system 

unit vectors in (§-, 7-, 6-) directions respectively 


angular frequency 
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Nondimensional quantities 
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a Us -/e 
subscripts 


(x, y) - component in (x-, y-) direction 


(§, 1) - component in (§-, T-) direction 
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I. INTRODUCTION 


Present generation computers make numerical approach to unsteady, 
two-dimensional subsonic flow attractive. In this study, two types of 
flow about a flat plate in nonlinear motion are investigated. Although 
having different objectives, related studies have been carried out by 
J. P. Giesing [1], N. D. Ham [2], and F. H. Abernathy and R. E. Kronauer ae 

Observations show that unsteady flow about a lifting body is 
accompanied by a vortex sheet being shed from the trailing edge. In the 
present study of unsteady flow about a flat plate, the plate is also re- 
placed by a bound vortex sheet, Both the bound and the free vortex sheet 
are then replaced by a finite number of discrete point vortices; thus, 
the integral equation of potential theory representing the field velocity 
is reduced to a finite difference equation. With appropriate application, 
this finite difference equation is used to determine the positions of the 
shed free vortices as a function of time and to express the boundary 
condition at any instant of time. In regard to the boundary condition, 
it is stipulated that the strengths of the vortices, both bound and shed, 
are such that the trailing streamline is tangent to the trailing edge. 

This satisfies the Kutta condition. Together with Kelvin's theoren, 
it is sufficient to determine uniquely the strengths of all vortices at 
any given instant of time. Steady flow is treated as a degenerate case 


of the above. 





The kinematics of the shed vortex sheet and strength distri- 
bution of the airfoil and free vortices are determined. Calculation of 
such important physical parameters as stability derivatives follows 


immediately. 


Il. THE BASIC EQUATIONS 


A potential flow model is used for unsteady flow about an 
infinite flat plate, Laplace's equation governs the velocity field 
and pressure can be calculated from Bernoulli's equation. From Kelvin's 
theorem the change of circulation about the flat plate is balanced by 
the shed vortex sheet(s)}]. Finally, it is assumed that the Kutta condi- 
tion is satisfied; this condition is sufficient to determine the rate 
of vortex sheet shedding. 

Using thin airfoil theory, the flat plate is replaced by a 
vortex sheet where vortex strength distribution is adjusted to make the 
flat plate a streamline of the flow. 

Consider an arbitrary flow field containing a vortex sheet 


A A 


with vortex strength distribution Y(s). Vorticity is considered positive 
A 


when inducing counterclockwise circulation. r is the position vector 
ed 





1. Bluff body flow is stipulated to be accompanied by vortex sheets 
shed from the leading and trailing edges. 





A 
on the sheet while s is the running variable along the sheet. The field 
A 
velocity at Ly is then given by 
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The same equation written nondimensionally is 
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Hereafter we will omit specific reference to the temporal argument. 
A description of the geometry of the field is appropriate. 


The absolute reference system is assigned coordinates (x,y,z) with their 


origin at the geometrical center of the flat plate. In the present study 


the center of the flat plate is therefore fixed but the plate is free 


to rotate in an arbitrary manner. i, j, and k are unit vectors in the 
Ped ~~ 


xX, y, and z directions respectively. A body fixed coordinate system 


(§, ll, S) with unit vectors e., e,, and e-, also has its origin at the 
oa | co 6 


geometrical center of the flat plate. 





Obviously, any point in the two-dimensional space considered 


may be described by either coordinate system. Suppose r is the vector 


from the origin to the point of interest. Then 
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bx + jy (3) 


Si oe cote (4) 


Furthermore, it is easily seen 
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Conversely, 
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If we decompose the 


and the free vortex 
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integral into the plate portion of bound vorticity 


sheet portion, we can write instead of (2) 


ds 
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A useful result obtained from (7) is the velocity normal to 


+ 
the plate at some point Be on the plate (7, = 0-). Thus 
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whereas the tangential velocity is 
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Now, for a flat plate whose angular position is given as a function of 
time, 6(T), the nondimensional boundary condition to be satisfied is 
Ve ( de { 1 
oe oS ae - --< € <-+¢e (10) 
dt Za ee 
8(T) is an arbitrary input to the calculation. By extending the boundary 
condition an indefinitely small distance downstream of the trailing edge, 


the Kutta condition is satisfied; i.e., the trailing streamline is tangent 


oe, 
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to the trailing edge. To determine the strength of shed vorticity, 


one makes use of Kelvin's theorem where 


v a 
—4{ /%s) dst = 0 (11) 
Dt 

The limit of integration can include any portion of the free vortex 
sheet in a Lagrangian sense or it can include the airfoil and any 


contiguous portion of the vortex sheet, 


TIT. THE FINITE DIFFERENCE SYSTEM OF EQUATIONS 


The numerical scheme to arrive at solutions for the airfoil 
and free vortex strength distribution and the kinematics of the free 
vortex sheet will now be discussed. As a reference, the flow chart (see Ap- 
pendix) will be instructive. The scheme is basically two-part. Knowing 


the strengths and positions of all vortices at t =t the velocity 


Re 
field is calculated and the new location of the free vortices are calcu- 
lated at t = th ook. The circulation strength of each free worees is 
conserved, Then the unknown strengths of the bound vortices and most 
recently shed free vortex are calculated by satisfying boundary conditions 
as described by the finite difference counterpart of (8) and (10). 

The vortex sheet representing the airfoil will be approximated 
by a finite number, N, of evenly distributed bound vortices. Bound 


vortex positions are indicated by circles in Fig. 1 (where, as an example, 


we choose N = 5), The distance between the vortices is 


A (12) 


O§ 





The matching points, eee at which the boundary conditions, equation (10), 
are satisfied will be halfway between consecutive vortices, so that as 

N- ©, the normal velocity will approach the principal value of the integral 
given by (8). Thus, the finite number of evenly distributed bound vortices 


will be placed at 


— sly . 
S = -O5 + var (so aecre N7 aes 


while the matching condition points will be located at 


L- O-5 

So = - O-5 + ret C2122). Ned) (14) 
The shed vortex sheet will be approximated as discrete doubly-infinite 
free point vortices shed from the trailing edge. A new vortex is shed 
at each increment in time and its strength is uniquely determined in 
satisfying the boundary conditions (equation 10). Once the strength of 
each vortex is determined, it will remain constant for all time. 

To start the calculation, let us suppose that the eer ear 
sheet distribution is known. The rate of change of position of the fe 


free vortex during a time increment At is given by the finite difference 


counterpart of (7). Replacing the vortex sheet strength YAs by ue at 


x,y = a yj 
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We have found it convenient to index by unity each free 
vortex, Thus the newly shed free vortex will always be identified by 
i = 1. Having established the positions of the free vortices according 
to (15) for t = th + At, combining (8) and (10) the boundary condition 
to be satisfied at that instant of time may be written in a finite 


difference sense as 
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with the right-hand side completely specified. (V is the component 


ie ip7 
of velocity normal to the surface at Sot induced by the fe free vortex. 


It is given by 


nr S 5; - (cos@ x. + SiNn@ y,) 
¥ | = -! ee : J , (18) 
digs 2k | (§. sine - ye + (§ . caso - Xi) 


It may be shown that Kelvin's theorem (11) may be written as 
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Equations (17), (18), and (19) may be combined and written as 
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b = TT (24) 


The coefficient matrix follows readily 


1 SE jose 2. ee 

ay = ee (25) 
2K Sj as Si ) = Le  N 

2 as | iv ny uy (t= 42,..,N) (26) 


Finally from (18) 
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Numerical representation of the boundary condition and Kelvin's theorem 
is now complete. Solution is obtained by inversion of the coefficient 
matrix. The equations were programmed in Fortran IV for the IBM 360-67 


(see Appendix). 
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IV. CALCULATED RESULTS 

The solution for the field about a flat plate in steady uniform 
flow is first considered. A comparison is made in Fig. (2) between the 
vortex strength distribution determined by the finite difference technique 
and that predicted by the exact integral (see, for example, reference baa 
Greater agreement is achieved by increasing N, the number of bound vortices 
defining the flat plate. This is, as noted previously, due to the fact 
that as N- ©, the solution from the finite difference technique approaches 
the principal value of (8). 

The case of impulsive acceleration from rest was considered 
mith o — 0.1 rad Geiser The resultant shed vortex sheet is depicted 
in Fig. (3) for T = 2.0 and tT = 3.0. ‘The vortex strength distribution on 
the plate approaches that of steady state asymptotically with time while 
the strength of successive shed vortices becomes negligible. Fig. (4) 
shows c, (normalized on angle of attack) vs. T. As expected, it approaches 
the analytical value of 27 asymptotically. 

Harmonic oscillation in pitch was then considered. 9(T) was 


defined such that 
a(t) = 9, sin WT 


The important physical parameters of 9, and W were set equal to 0.1 and 
0.1, 1.0, 2.0, and 3.0 respectively. The number of bound vortices 
defining the plate, N, was set to 11 and the maximum number of time steps 
taken was 150. AT was fixed so that the longitudinal spacing of shed 


vortices would be of the order of A&, 
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Figs. (5-10) show the resultant shed vortex configurations 
for the specified frequencies at selected elapsed times. Since the wake 
is approximated by discrete free vortices, curves drawn sequentially 
through the points yield an approximation to the shed vortex sheet. 
As distance downstream increases, definition of the sheet becomes more 
aiecficuit,. 

The short period frequency of a stable aircraft at cruise 
speed corresponds to a reduced frequency of the order of 0.1. As seen 
in Fig. (5), the resultant wake has negligible vertical deflection. 
This reaffirms the assumption that the wake is linear for most practical 
purposes in unsteady aerodynamics. Comparison of Figs. (6-10) demonstrates 
the increasing nonlinearity of the wake with increasing frequency. 

An interpretation is made in Fig. (10) of the vortex sheet 
for © = 3.0 and T = 10.0. At a sufficient distance downstream, the wake 
rolls up periodically into concentrations of free vortices. The spacing 
ratio (that is, the ratio of the vertical displacement of the "centers of 
gravity" of cloisters of vortices of opposite sign to the horizontal 
distance between cloisters of like sign) is typically 0.07. The classical 
ratio predicted by von Karman [5] for a stable configuration of double 
infinite rows of point vortices of alternating sign is 0.281. 

Coefficients of lift and drag were calculated by applying the 
Kutta-Joukowsky Law locally to each bound vortex defining the flat plate. 
Figs. (11, 12) present these values for W = 3,0 as a function of elapsed 


time T, A typical value for cy, at T = 0.50 and @ = Soil. is 22 ek 


eum —“ag> 
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The steady state value for this particular angle of attack is .6l. 
Cc, and cg are periodic with time with negligible phase shift relative 


to the harmonic oscillation. 


V. CONCLUSION AND DISCUSSION 


The method developed herein yields good results for the unsteady, 
two-dimensional flow about a flat plate executing two separate nonlinear 
motions: impulsive acceleration ard sivall amplitude oscillation about a 
Spanwise axis. With the flow field specified, such significant quantities 
as stability derivatives are obtained readily. 

In future efforts, bluff body flow, such as large amplitude 
oscillation or tumbling, with vortex sheet shedding from the trailing and 
leading edges, is contemplated. Preliminary investigations were made. 

Once again the potential flow integral equation representing the field 
velocity (1) is replaced by a finite difference counterpart. The boundary 
condition, (8) and (10), is extended an indefinitely small distance upstream 
of the leading edge and downstream of the trailing edge. (15, 17, 18, 19, 20) 
are modified accordingly. Kelvin's theorem determines the rate of vortex 
shedding and the two Kutta conditions establish uniqueness, Initial attempts 
to establish the flow field were unsuccessful, These attempts suggested 

that further refinement is needed in the numerical model. Consideration 

must be given vortices in close proximity (possibly in the form of the 
criterion proposed by Ham [2]). Furthermore, some modification may be 


required for the boundary condition at the extremities of the lifting surface. 
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APPENDIX 


OSCILLATING FLAT PLATE 


COMPUTER PROGRAM 





ALCHS Sa, 


UTAIS TAIT OFPTALIZO NO 


MARZQOFL Arla) 


SPECIFY INITIAL CONDITIONS , 


CALCULATE STEADY STATE FIELD 


SPECIFY U_(T), 6(r) 





CALCULATE VELOCITY FIELD 


Cet = t+ Ot | 


CALCULATE POSITIONS OF FREE VORTICES 


CALCULATE STRENGTHS OF BOUND VORTEX SHEET 
AND FREE VORTEX MOST RECENTLY SHED BY 
SATISFYING BOUNDARY CONDITIONS 


PRINT 





MOO 


AOD OY CHO ONMOKC KONO 


OOOO 


Pe oe ae 


L3G 


140 


PUNCH D&eCK, OSCILLATING FLAT PLATE, 


REAL 4NKXM NYE 

CpMiews iON EtLil, ll}, Ft ll), NGI11), NPULL), Wtil), CULL) 
DIMENSION ABC L1), YBC(Lils XF( 200), YRF(Z00), AMCIL1), YNC(11) 
CiMENSION GAMMAF( 200), ALM(11) 

DIMENSION All2,12), B12), Cll2), L&UL2), NVI L2)d, VNE(1C) 
DIMENSION TGAM( 200), ZX(200), LZY(200), X18B(11L) 

Pl = 6228318 


fees THE eNUMBER OF EVENLY DISTRIBUTED. BOUND VORTICES CN THE LIFT— 
MiworsoURPFACE. “AN 1S THi INCREMENTAL DISTANCE TN The FINITE OCLFFE- 
eeNvCc SCHEME. 


T IS THe ELAPSED TIMt, TV TS THe TWNITIAL TIME AND DELT 15 TRE Tie 
INCREMENT. 


Poetic lS THE AMPLITUDE OF GSCILLCATICON IN RADIANS. WHILE FRE IS 
The REDUCED FREQUENCY OF OSCILLATION. 


Roa 1 

Aii = 1.0/FLUAT(N-1) 
SiS N=1 

AeA NE 1 

GHetAZ = 0.1 

FREQ = Qel 

ier. OO. 1 

eee ene 


TZ =7- 71 


Sar 307.5 = Ly 

ALB(J) = -0.5 + AN * FLUAT(J—-1) 

AIM(J) = -0.5 + AN * (0.5 + FLOAT(J—1)) 
CONTINUE 


THETA TS Tine MGUIFIED FREE STREAM ANGLE WHILE CineTA [IS THE CHANGE 
PNetnelA wihih KResrpec hh aG TlME. 


TRETA = THETAZ * SIN(FREG * TZ) 
DeGREE = THETA * 57.29583 
DIKETA = THETAZ * FRLG * COST FREQ * TZ) 


COSTi = COS( THETA) 

SINTH = SIN(THETA) 

Ale = OwS * COSTH 

YTc = 0.5 * SINTH 

XXX = SCRTCATE * XTE + YTE * YTE) 
NAM = -YTC/SXKX 

NYM = XTESAXX 


DE 0 Ce 0 es ae Ger 


Rots) = ALStI) * CUSTi 

YoCJ) = X1B(J) * SINTH 

AM(J) = X1M(J) * COSTH 

YM(J) = KIM(J) * SINTH 

CUNTINUE 

CO 150 [IT = 1,N 

GOO 150 J = LyN 

Eliya) Sa ee? Wee TOL A Gal) A=" Oe 90.) 


CONTINUE 





a ak ils i al a 


B50 


DG 3907] = "PN 

FCI) = NAM 

CONTINUE 

IXRUw = 1 

CALL MINV(EsN,gDO,NGyNP) 
CALL GMPROCE,Fy,HsNygN,TRUW) 
Wee GO Le<= Ll» W 

oe.) = Hi 1) 


60 CUNTINUE 


MOON al al a ia MOO iad 


ll 


OOM 


2V0 


mein it ITAL BOUND VORTEX STRENGTH DISTRIBUTICN, 13 SPCCIFIEG. 


CONTINUE 

PONS TEADY FINITE DIFFERENCE. SOLUTILN. FULLOGKS. 

Boe LOCO M=1, 150 

je)! + DELI 

TZ = T - TT 

Mme eReEt VORTICES ARt ASSUMED SHED FPROMFTHE TRAILING ELGE-. 


XFC 1) 
Were Gl} 


ABC) 
YEN) 


TRE VtELUCITY FIELD IS BeETeRMINED. 
CGO 800 J = ly M 


WerOCily INDUCED BY THe [Th FREE VORTEX GRIGINSLLY SHED FRCY Ere 
MroiLinNGeeeGe ON THe JTH FREE VGRTIEX GRIGINALLY SaeD FREM [at 
TRAILING EDGE. 


wot = C.0 
VYTT = ULC 
IPt“.£0.i) GO TO 5O0l 


GAMMAF( 1) IS AS YeT UNKNOWN. 
DG 500 I = 2,ytt 
A VORTEX CANNOT ACT ON ITSELF. 


IFCT-&E@.4) GO TO 500 

Keeani sho AECL) 

Mao Nr = YE) | 

KA = K ¥* K 

Vie = yy oY 

Ewin: == Pla AK + YY) 

VATT ((-Y¥ * GAMMAF(CT))/DENCM) + VXTT 
VYTT (( X * GAMMAF(1))/DENGHM) + VYTT 
CONTINVGE 


ito ot 


DOS COUN Tse 


mer CC ans 





_— 


Or etrerCT oy OO.) 


oy CTO) 


So 


= ty — 


600 


~_ 


CY CPO 


| 600 


Ire). CMC Ces 


860 


aay 


VELOCITY INDUCED BY THE ITH BOUND VORTEX ON THE JTH FREc VORTEX 
URIGINALLY SHED FROM THte TRAILING tOGE. 


VAET 


C.Q 
VYbBT U 


O 


Won 


THE BOUND VURTEXK AT XB(N) ITS EACLUDED TO ARTIFICIALLY SUPPRESS 
THE SINGULARITY. : 


CO600 1° = Ly, 

Met let wsNeANDsedJsEQs<1) GU 16° 600 
KX = XF(J) —- ABIL) 

Meee YF(J) -— Y6(1) 


KK = K * X 
iore= Y ~ Y 
wenn = Pl *. (XA + YY) 


eat = ((=¥ + BLE} * AN) /DENGM) + VABT 
VY¥8T = (CC X * Bl) * AN) /DENCM) + VYBT 
CONTINUE 


Pome eoli ino OF THE FREE VORTICES ARE DETEXRMINeEL (PLN 
ie PUK) oF OODELT. 


ma) = AVALIN+*VABT) * DELT + DELL + XE 62) 
Z¥(J) = ( VYTT + VYBT) * DELT + YF(J) 
CONTINUE 

Ou 300 J = i,i 

XF(J) = LX) 

Vitak =. 2Z¥ (3) 

COGNTi NUE 

ThetA = Thevad * SINCFREQ * 72) 
DEGREE = THETA * 37-29983 

ALPHA = -—(DEGREE) 

CTHETA = THETAZ * FREQ * COS(FREQ * 12) 
COSTH = COS( THETA) 

SINTH = SINCTHETA) 

ATE = 0.5 * COSTh 

Yit = 6.5 * SINTH 

KEG = -~ATE 

Yiee, = —YTE 

KAN = SOGRVCATE * XTE 4+ YTe * YiE} 

wa = =YTEG/XAAX 

NWYM = XTE/SXXX 

DU 230 J = 1,N 

ROG Aol) * CUSdk 

Yu(J) = xXiBl(J) * SINTH 

AM(J) = X1IM(J) * CUOSTH 

YM(J) = ALM( J) * SINTH 

CUNTINUE 





SR a a a i a a a ol a a a 


al aad Si al 


1Q0 
OE 


200 


i il Sa taal Sa 


620 


CVO OOOO 


LLO 


Nas 


al a 


meeCcUuLATING THE COMPONENTS UF THe ROW VECTUn, Cll). 


BACH ELEMENT OF THE RUm% VECTOR ITS, AT A PARTICULAR FMAICHING 
Combi TiN POINT, AWN EXPRESSIUN GF THE NCRMAL CCMPUNENIS UP FREE 
meee VELOCRIY AWD VELOCITY 1NDUCCD BY FREE VERTICIES (CEACEPT 
Meesc MUST ReCENTLY SHED) PLUS WHAT COMPCNENT CUE TC RCOTATICHN CF 
Wee LIFTING SURFACE. 


Wes Z200 I = ign 
PAECULATING VHE INCUCED VELOCITY AT THE ITH MATICHING CCHNOITICN 


POINT OUE TO THE JTH FREE VORTEX ORIGINALLY SHED FROM THtE TRAILING 
Buccs Wilh lie EXCEPEIUGN OF THAT UNE MOST RECENTLY SHEL. 


VX = CeO 

wie O60 

teeitiscvel) GO 10° 101 

UG 100 J = 2, i 

K = XM(1) — XF(J) 

feet Gl) YP) 

AK = K * X 

Vor oF o¥ Y 

eNom “= —~P1 * (AK + YY) 

VA = ((0-Y * GAMMAF(J))/0ENCM) + VX 
VY = (CC X * GAMMAF(J))/0ENOQM) + VY 
CUNTINVE 

CONTINUE 

VinM( TD) = (VX * NAM + VY * NYM) 


Cll) = NXM # VNM(I) -— (xXIMCL) * DTHETA) 
CONTINUE | 


CermPONENT DUS To KELVIN*S THEOREM. 


SUrGAM = 0.0 

UU o20 I = 1,N 

SUMGAM = o(i) * AN + SUMGAM 
CONTINUE 


C(N2) = SUMGAM 
CONSTRUCTING THE COEFFICIENT MATRIX. 


THOSE ELECMENTS DUE 10 THE VELOCITY INDUCED BY THe JRO BGUNG VORTEA 
AT tHe JTR MATCHING CONDITION POINT. 


We 2G Ae = laa 

Cu LlC J = ly WN 

Migkgo) = Ale/(r lee CELGAT( I-11) = 0245)) 
CONTINUE 

CONTINUE 





Stotetatamiaet 1.0 ks oC 


Pee eee timo DUE STO IHE crreGl OF THE PREE VORTEX MUST RECENTLY 
SHED FROM THE TRAILING EDGE CN THE ITH MATCHING CGROITICN PCINT. 


Seo = Pl * (Xx + YY) 
; = -~Y/DENCM 
| Vy = X/CDENOM 
ACIyN+1) = -(NXM * VK + NYM * VY) 
530 CONTINUE 


C 

é Mmeoe CLEMGH1S DUC 10 COUNSTDERAVIGN OF KELVIN'S THEGRor. 
& 

| 


CU 550 J =1,N 
A(N2,J) = AWN 
550 CONTINUE 
| A(N2,N2) = iL. 


C PNVER STON GF Al1»yJ?) « 
iC Pivewou (Adi yJI)) * Cli) = b(d) 


SCAU GMeRD(A,C,6;N2,N2,1R0W) 
GAMMAF(1) = 3(N2) 
[Pte eC.150) UD 10 9S9 
Geir) G59 
999 CGNTIWNUE 
—wWkITE(7,20) Ty FREQ, ALPHA 
20 FORMAT(2r 5.2, FLO.5) 
WRITE( 7,29) XLEy YLE 
25 FURMAT(2F10.5) 
WRITE( 7,260) XTE,y YTE 
26 FORMAT(2F10.5) 
Mollet 7,21) (XE 0)y YE (Jey Jal) 
21 @PGAMAT(2F1i0.5) 
WKITTE(7,22) 
22 FERMAT(C!1COG.") 
699 CUnTINUE 
TGAM(C 1) = GAMMAF(I) 
Pot) ae 
Dy) = Vet) 
TCO CUNTINUE 
DO 7l0 1 = 1,)h 
GAMMAF(I+#1) = TGAM(T) 
APL Ly = «£ ) 
Yetiel) = Z ) 
710 CONTINUE 
L0CG CONTINUc 
STUP 
END 
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